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ABSTRACT 


A technique is developed to determine the rotational 
temperature of Nitrogen Molecular Ion (N^) from the emission 
spectra of B-X transition, when P and R Branches are not 
resolved. Its validity is tested on simulated spectra of the 
0-1 band of N 2 + produced under low resolution. The method is 
applied to experimental spectra of N 2 + taken in the shock 
layer of a blunt body at distances or 1.91 cm, 2.54 cm, and 
3.18 cm from the body. 

The laser induced fluorescence (LIF) spectra of copper 
atoms is analyzed to obtain the free stream velocities and 
temperatures. The only broadening mechanism considered is 
Doppler broadening. The temperatures are obtained by manual 
curve fitting; and the results are compared with least square 
fits. The agreement on the average is within 10%. 
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INTRODUCTION 


There is a strong need for non-intrusive diagnostic 
techniques such as passive radiation techniques which utilize 
radiation emitted by the plasma or active radiation 
techniques such as Laser Induced fluorescence (LIF) 
measurements. These techniques are essential to the arc get 
flow studies in order to fully understand the nonequilibrium 
flow conditions spacecraft encounter during reentry. The arc 
jet facilities simulate the spacecraft reentry conditions and 
are used for testing thermal protection materials. The 
identification, characterization, velocity and temperature 
determination of different atomic, molecular and ionic 
species makes an important contribution to the arc jet flow 
studies. The major species of interest in the spacecraft 
reentry flow environment are N 2 , 0 2 , NO, N 2 + , N, O, 0 2 , N0 + , 
N , 0 + and electrons. Concentration of these species change 
when the enthalpy of the flow changes. This report is divided 
into two parts, one on the analyses of the emission spectra 
and the other on the analysis of the LIF data. 

Emission studies: During the last couple of years the arc 

jet flow diagnostic program at NASA, Johnson Space Center 
reported emission studies (Ref.l and 2) of N2 and N 2 for 
vibrational and rotational temperature determination across 
the shock layer. This technique involves calculating the 
spectrum for a number of cases and obtain integrals over the 
wavelength regions of the spectrum as a function of 
temperature. Ratios of these integrals were then related to 
the temperatures used to generate the spectra. Spectral 
integrals from the measured spectra are then compared with 
the calculated values to determine the temperature. This 
technique has some limitations and requires numerous 
parameters to produce the calculated spectra. This method 
requires repetitive calculations be performed each time a 
comparison is made. In this report a simpler technique is 
presented to find the rotational temperatures of N 2 in the 
shock layer. This can be applied to other molecules as well. 

LIF studies: The LIF measurements reported here probe 
different species with a known excitation mechanism using 
selective excitation with a narrow band laser. Preliminary 
LIF studies of copper are used to calculate the velocity and 
temperature are presented here. 

MEASUREMENTS 

Emission spectra: Details of the experimental set up to 
record the emission spectrum of nitrogen in the wavelength 
region of 340 run to 480 nm are reported earlier (Ref. 3). 
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LIF studies on copper: Laser fluorescence spectra of copper 
recorded earlier using the experimental setup shown in Fig.l, 
are analyzed. The setup essentially consists of a YAG pumped 
dye laser system with intracavity etalons in the dye as well 
as in the pump laser reducing the dye fundamental band width 
to less than 0.08 cm -1 . To probe the copper atoms, 
Dicynomethylene (DCM) dye is used together with a KD*P 
doubling crystal to produce laser light at 327.5 nm. Details 
of the laser fluorescence detection and control electronics 
along with laser pulse energy used for copper are reported 
earlier (Ref. 4). Velocity measurements of the free stream for 
various mass flow rates and arc currents have been carried 
out using copper as the tracer material. The excitation 
wavelength is about 327.5 nm (Figs. 2) and the 578 nm 
fluorescence is collected at right angles to the flow using 
cutoff and narrow band interference filters. Two separate 
measurements were taken for Doppler shift studies one with 
the laser beam perpendicular and the other at an angle of 
60° with respect to the flow direction. 

RESULTS AND ANALYSIS 

Emission spectra of N 2 + : The rotational structure of the 0- 
1 band of B-X transition of N 2 + is shown in Fig. 3. Due to 
the limited resolution the P and R branches are not resolved. 
In order to analyze and determine the rotational temperature, 
we have developed a method described in an earlier report 
(Ref. 3). For completeness of this report the equations are 
shown in Appendix 1 and 2. In order to test the accuracy of 
these equations, simulated spectra of N 2 using NEQAIR 
program under low resolution are produced at various 
temperatures. A sample spectra produced at 8000 K is shown in 
Fig. 4. We have now used the areas and corrected them for P 
and R branches using equations (7) and (8) shown in Appendix 
2., at various temperatures (T) . Graphs are now plotted: In 
[Ipc/ (k'+k"+l) ] against k* (k'+l) for P-branch and In r 
Ir c / (k'+k"+l) ] against k' (k'+l) for R-branch. According to 
equation (4) shown in Appendix 1, as given by Herzberg 
(Ref. 5), the slope of the line gives B 0 hC/KT from which the 
rotational temperature can be determined. The temperatures 
(Tp and Tp) are now calculated using the slopes of the 
straight lines obtained from the graphs for both P and R 
branches. The particular set for which the temperature Tp = 

T R = T is considered as the rotational temperature. Fig. 5 is 
a graph showing the agreement between P and R branches at 
8000 K. Table 1 shows the temperatures calculated using the 
slopes of the lines at three different temperatures. 

This method is now applied to the spectra taken on 
February ' 89 at 1.91 cm, 2.54 cm and 3.18 cm positions in 
the shock layer of a blunt body . Figs. 6 and 7 shows the 
graphs plotted for corrected P and R branch intensities 
taking the areas and the results are tabulated in Table - 2. 
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The calculated temperatures are slightly less than the 
temperatures obtained by Blackwell etai. (Ref. 4) using 
computer simulation techniques. This low temperature seems to 
be due to the overlap of intensity of some other bands of 
the same species or due to the emissions of the species 
present in the shock layer in the region analyzed. If an 
account can be made for this background intensity rotational 
temperatures can be more precisely determined even in 
unresolved bands. Work in this direction will be continued in 
the future. 

LIF studies on copper: A representative spectrum of Doppler 
shifted copper spectrum is shown in Fig. 8. The measured 
Doppler profiles are used to determine the velocity and 
temperature. The velocity is obtained using the relation 

V = CAW/ W Q (1) 

where V is the velocity of the free stream, C is the speed of 
light , W is the measured Doppler shift and W Q is the laser 
excitation frequency. A graph of the measured velocity 
against bulk enthalpy is shown in Fig. 9 and it demonstrates 
the expected increase in velocity with increasing enthalpy. 


To determine the temperature, a linear curve fit of the 
quadruplet of Gaussian profiles was fit to the measured 
profiles. The laser line width was accounted for by folding 
its profile with the Doppler profile. Both were taken to be 
Gaussian. 

The intensity of Gaussian line shape is given by 

I (W) = I 0 exp[ -4 ln2 { W Q - W/AW d > 2 ] (2) 

and the Doppler Width & W D is given by 

&W D = w o/ c ( 8 RT ln2/M) V 2 

where 

I c = the peak intensity of the Doppler profile 
W Q = the frequency at peak intensity 
R = gas constant 
T = Temperature in Kelvin 
M = Mass of copper atom. 

Fitted curves at different temperatures are developed 
using equation (2) and were overlapped on the experimental 
curve. The one that best fits to the experimental curve is 
taken to be the temperature at that flow rate. A sample curve 
fit using one of the two peaks is shown in Fig. 10. The first 
peak is saturated and so only one peak is considered for an 
accurate determination of the temperature. Temperatures are 
determined using this method at various flow rates and the 
results are presented in Table - 3. The results do not 
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compare particularly well with the least square fit, 
indicating the assumptions used in the simple calculations 
are not adequate. More systematic study has to be carried 
out and using least square fit Gaussian profiles may yield a 
more precise temperature determination. It may be noted that 
a Voigt profile (a combination of Gaussian and Lorentzian 
profile) may be the appropriate choice for future studies. 


CONCLUSIONS 

On the basis of the analysis made in both the projects the 

following conclusions are made: 

1. Rotational temperature measurements can be made within 5% 
even in the case of unresolved bands due to a single 
species if there is no overlap of other bands. 

2 . This method can be extended to the observed experimental 
spectra, where other species are present provided the 
contribution due to the background is accounted for. 

3. The method used to find the temperature of Doppler shifted 
curves in copper agrees on the average with in 10% with 
that obtained using least square fit. The velocity can be 
calculated to an accuracy of 10%. 
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TABLE 1 


SIMULATION SPECTRA - COMPARISON 


T 

(K) 

(!> 

% 

Error 


(S) 


% 

Error 


8000 

7906 

1.1 


7553 


00 

• 

in 


6000 

6234 

3.9 


6388 


6.5 


3000 

3125 

4.2 


3217 


7.2 


TABLE 2 

.- ROTATIONAL TEMPERATURES 

OF 

SHOCK 

LAYER 

Position of 
Shock layer 
(cm) 

% 

% 

Error 

(Jb 


% 

Error 

1.91 


3991 

0.4 


3982 


0.2 

2.54 


2954 

1.5 


3045 


1.5 

3.18 


6337 

5.0 


6131 


2.2 



TABLE 3 

.- COPPER 

DOPPLER 1 

SHIFT DATA 


Current 

(Amps) 

Flow Rate 
(lb/s) 

Doppler 

Shift 

10” 3 nm 

Velocity Temp 

(m/s) (K) 

Curve least 
fit Square 
fit 

Enthalpy 

(Btu/lb) 

700 

0.05 

1.486 

2721 


»» 

7000 


0.10 

1.500 

2747 

— 

2022 

5640 


0.15 

1.480 

2711 

1200 

1148 

4900 


0.20 

1.468 

2688 

1500 

1486 

4150 


0.25 

1.519 

2782 

— — 

1076 

3760 

500 

0.10 

1.547 

2834 

1400 

1457 

4260 


0.20 

0.885 

1621 

1300 

1478 

3100 

300 

0.10 

1.232 

2256 

1100 

1117 

2780 


0.20 

0.709 

1299 

600 

498 

1925 
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APPENDIX -1 


ROTATIONAL TEMPERATURE EQUATIONS 


I(P - Branch ) = A - 2 Je 


B 0 J(J - l)hc 
kT 


( 1 ) 


I(R - branch ) = A - 2 (J + 1 )e 


B 0 U + 1)(J+ 2)hc 
kT 


( 2 ) 


Rearranging equations (1) and (2) 


In 


2J 


. i. * - !• ” - 1 »“ 

kT 


In 


Ip 
2 JJ 


B 0 J( J — l)hc 

kT 


(3) 


Similarly for R Branch 


Ini 


l 2 0+ DJ 


B ft (J + 1)(J+ 2)hc 

_ A • o 

' A ' kT 


(4) 


Equation (3) and (4) can general be written for 
either case as 


1 J I ' 1_ A Bl'(k'fl)hc 

in L k' + k" + lj" A " kT 


where k‘ - upper rotational state quantum number 
k" - lower rotational state quantum number 
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APPENDIX - 2 


CORRECTE D, ..IMIEMS II IES 


*exp *p + *r 


( 5 ) 


EXP 


+ 1 


( 6 ) 


Let 


I 


B where L and are the 


intensities calculated for any temperature 
using equations (1) and (2) 


* Rc * EXP 


L1 + BJ 


( 7 ) 


*Pc ~ ^EXP 


1 


1+ i 


( 8 ) 


I Rc and I pc are the intensity contributions 
of P and R branches in the observed spectra 



